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Abstract 

We study the contributions of the CP violating anomalous WW^ interac- 
tions to 6 ^ sl'^l~ . We obtain cutoff independent results on k and A, by 
constructing an asymmetry for the process B — > K*£~^i~, where the B and 
B events are added. We show that a sample of 10^ B K*e~^e~ events can 
yield a bound, \k\ < 0.42 at 90% C.L., which is much tighter than the recent 
constraint from DO. 
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Two of the aspects of the Standard Model(SM) which remain to be further explored are 
the gauge nature of the W, Z bosons and CP violation. The vector boson self interactions 
are uniquely fixed by the SU{2)l x U{1) gauge structure of the SM. However, new physics 
as well as higher order corrections will modify these self interactions. While tests of trilinear 
gauge couplings are now being done at present colliders, only weak bounds on them exist 
imj^. CP violation on the other hand is only parameterized in the SM via the CKM matrix, 
and has been seen only in the neutral K mesons. Much effort is now being put into studying 
CP violation in the B mesons. Large number of B mesons are expected to be produced 
in the near future, at the new asymmetric colliders and upgraded Tevatron, which would 
enable study of its rare decay modes. Rare decays of B mesons, occur in the SM only through 
loops and are thus expected to be very sensitive to new physics. In particular the mode 
h st^i~ has been extensively studied for possible contributions from new physics. 
New physics that affects the gauge nature of the W, Z bosons, can be parameterized in terms 
of CP conserving as well as CP violating anomalous trilinear gauge boson couplings. The 
CP conserving couplings have been considered in the rare modes h — > S7, h — ^ sl'^fr l^-Q. 
Since, the origin of CP violation is not understood, any potential source of CP violation 
should be pursued. Here we focus attention on the CP violating anomalous couplings, by 
studying the mode B — > K*i^i~ . 

The couplings of W bosons to the photon and Z can be described by an effective 
Lagrangian, the most general form of which may be written down with the minimum 
requirements- that of Lorentz invariance, global SU{2) and local U{1) symmetry |^ as: 
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jC'^ff = - 9v [^9Y {wipW' - W^'^W^p) + iKyW^WpV 

In the above equation V represents the neutral gauge bosons either the photon or the Z, 

V^p = d^Vp - d/sVa, W^p = d^Wp - dpWa, V^"^ = ^e'^i^P'^Vp^ and gv is the WWV coupling 
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strength in the SM with g-y = e and gz = ec/s, where = 1 — = M^^/M'^- In the SM, the 
couphngs gi, Ky are unity and all others are zero. New physics may result in a modification 
of these couplings. The deviation from the SM, the so called anomalous couplings, need to 
be constrained experimentally. We shall here concentrate on the CP violating couplings k 
and A. While many theoretical and experimental studies of anomalous triple gauge boson 
couplings have been done, few bounds exist on the anomalous CP violating couplings. 
Direct bounds have recently been obtained |1| on k and A from pp WjX at the tevatron. 
Stringent indirect constraints [|TU[ on these couplings come from electric dipole moment 



(EDM) of the neutron and electron. However, these constraints assume naturalness and 
are cutoff dependent. Hence, bounds from CP violating asymmetries that we pursue here, 
would be complementary to the EDM bounds. Possible contributions of the couplings k and 
A to processes at colliders have also been examined for LEPII |^, upgraded tevatron |Tl 



as well as future linear colliders []T2|. In the chiral Lagrangian approach the dimension 6 
operator with coefficient. A, is ignored. For '6' penguin modes such as 6 ^ 57 or 6 ^ si'^i~ , 
unitarity of the CKM, results in cutoff independence for contributions from k and A. Since 
both k and A contributions are finite, requiring no cutoff, we do not neglect A. However, 
the constraints on A are expected to be weaker. 

The effective short distance Hamiltonian relevant to the decay b s£^£^ |13- |T5|] leads 
to the QCD corrected matrix element, 

V27r 

(2) 

where only the dominant top quark contribution to the loop is retained. Cj {j = 7, 8, 9) 
are the Wilson coefficients given in Ref. [l^JT^, is the mass of the b quark, is the 
invariant lepton mass squared, b^ ji = (1 =i=75)/2 6 and Vt = V^^Vtb is the product of the CKM 
matrix elements. The anomalous WWV couplings, result in a shift in the values of the short 
distance coefficients C/s at = Mw, 

Cj = C/*' + I kC^ + I ~XCf, (3) 



where, C^^ are the SM coefficients and Cj^, C'^ are the contributions from the anomalous 
couphngs. These shifted values of Cj's are then evolved down to the h quark scale. The 
anomalous WW-^ coupling coefficients have been evaluated to be, 



_ ^\ _ ^ ( 2 + 37 Xt + lOx^ - Xt (3 + 5xt) \ nxt^ 
^" ^"Ma^l2l, 6(l-a:.r + {1 - Xt]' 



n'K ( 1 ^ Xt {3-Xt) \nxt \ , . 



^ _xt [l- xf + 2xt \nxt) 

— ~ " 



2(1 -Xtf 



where, Xt = y7T^ mass of the top quark. The WWZ coefficients can be related to 

the corresponding WW''y coefficients listed above, by modifying the couplings and propaga- 
tor. Cg'^ also contributes to the WWZ vertex, with a functional form in xt, identical to that 
of Cg'^. However, the anomalous contributions from the WWZ vertex may be neglected in 



comparison to that of the photon. Our results for are in agreement with those of refs 



and [H] but we disagree with ref. M on Cj and on Cg . Cg has been evaluated for the 



first time. Note that Cg turns out to be exactly equal to Cg , perhaps indicating that for 
the contact term Cg, the derivatives from dimension 6 operator and the gauge part of the 
propagators contributing in case of the dimension 4 operator, behave similarly. Another 
interesting feature is that Cg is doubly chirally suppressed for the CP violating couplings 
k, A. 

Previous studies |^J^ of CP violating anomalous gauge boson couplings in rare B decays, 
have used the mode b — > 57, to constrain k and A separately, disallowing any possible can- 
celations between them. They also used only the decay rate of 6 — > 57 and hence depended 
quadratically on these couplings. In order to achieve maximum sensitivity to CP violating 
couplings and to have a clear signal of CP violation, we need to estimate observables that 
are explicitly CP odd and must therefore depend linearly on the CP violating anomalous 
couplings. CP violating asymmetries have also been considered in 6 57, in SM and two 



Higgs doublet model |T^-|T8[. However, such asymmetries rely on the presence of a large 
strong phase arising out of final state interactions. CP violating asymmetries in rare modes 
that require flavour tagging are going to be very difficult to detect. We therefore prefer 
to use a technique that neither needs flavour nor time tagging. Such an asymmetry has 
been considered in ref. We adopt this technique and construct asymmetries for the 

mode B K*i^i~ non resonant), to obtain simultaneous bounds on k and A. The 

mode B K*i'^i~ is theoretically clean and can provide additional information due to its 
richer kinematics. We present here only the case of charged -B's. For neutral -B's, after time 
integration, the asymmetries that we construct, have exactly the same form as that for the 
charged B's 

The transition matrix element for the exclusive process B{p) —>■ K*{k)£~^£~ — >■ 
K{ki)Tc{k2)i^{qi)i^{q2) can be written for each of the operators in Eq.(||) as, 

{K7i\Ma^,{l ± i5)q''b\B) = zAe^^^^K'^k'^q^ ± BK^ ± Ck^ , 

{K'k\s^^{1 t l5)b\B) = tVe^^afsK'^k'^q^ ± SK^ ± J^k^ . (5) 

The form-factors ■ ■ ■ , JF are unknown functions of q"^ = {p — k)"^ and other dot products 
involving momentum, k = ki + k2 and K = ki — k2 and can be related to those used in , 
as given in the Table. The variable cr, in the table, arises due to the decay of K* Kn, 
evaluated in the zero width approximation. They can be similarly related to those in heavy 
quark effective theory (HQET) The HQET form factors cannot currently be reliably 
predicted over the entire dilepton mass range. The results that we shall obtain, do depend 
on the numerical values of the form factors. In future, with large number of -B's available, 
the form factors can be determined experimentally. The current proportional to does 
not contribute as it couples to light leptons. In our notation MB,'mK*,mK and are 
the masses of the B, K*, K mesons and the pion respectively. The matrix element for the 
process B — >• K*i^i^ — >• Kni^i' can be written as 

M{B ^ Kntr) = '^{(iaLe^,^pK-k''q^ +f3L + Pl k^) L £ + L ^ r} , (6) 
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(1 ip ■yA 

where L,R = , q = qi + q2 and Q = qi — q2, and the coefficients ai^ji, [3l^r and pl^r 



are given by 
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(7) 



In the above equation a, (3 and p are recast in terms of a, h and r so as to identify the strong 
phases 5 and the weak phases 0. Using CPT invariance, the matrix element for the decay 
B Kfri^i^ can be obtained from the B Kni'^i^ by replacing aL,R —CiL,R, I3l,r 
pL,R, PL,R Pl,r HES' where 



aR,L = exp {i^'^^i) exp (-^0^,^) (8) 

and similar relations hold for (3 and p. The matrix element mod. squared for the process 
B K*i'^i'' — >• Kni'^i" is worked out retaining the imaginary parts in a, [3 and p to be, 



\M{B{p) -^K{h)7r{h)i+{qi)r{q2)) 



2 a^Gl 



27r2 

{ 2 t.^p^k^^K'^qPQ" {K ■ Q lm{aL {31 + an {3*^) + Im(p^^ {3*^ - p^ PI) 
-k-Q Im(pi al + p/j a^)) + 2 Re(pL al - pR a^) 
(-k ■ K q^ k ■ Q + k ■ qk ■ Q K ■ q + m\, q^ K ■ Q - k ■ q^ K ■ Q) 
+ 2Re(pL/?2 + pR/?^) [-k ■ Kq^ + k ■ qK ■ q - k ■ Q K ■ Q) 
+ 2Re{aRpR-aLPl) (k^ q^ k ■ Q - k ■ Q K ■ q^ - k ■ K q^ K ■ Q 
+ k-qK-qK-Q) + (pl + pl) Um%, q^ + k ■ q^ - k ■ Q^' 



(9) 



where 



+ (ai^ + aR^) (-K^q^k-Q^ + k-Q^K-q^ + 2k-Kq^k-QK-Q 

-2k-qk-QK -qK -Q- q^ K ■ Q"^ + k ■ q"^ K ■ Q^' 
+ {PL'+PR'){-K'q' + K.q'-K.Q')} 

k ■ K = m% - ml, q - Q = 



6 



k-q=]^{Ml-ml.-q^) 
k-Q = XMBCOsOe 

2 2 

K-q = \K XMb cos Ok* + k ■ q 



K ■ Q = k ■ Q 7^ + Xk (k ■ q cos 9i cos 6k — \/ q'^m\, sin 9i sin 9k cos ^p) . 

Here X is the three momentum of the i^i^ or Kit invariant system in the B meson rest 
frame and is given by 

^ ^ {k ■ q^ - q^m]^,)l 
Mb 

Xk is related to the K three momentum in the K* rest frame and is defined as 

and similarly, Ae is related to the lepton momentum in the i^i^ rest frame and is given by, 



g2 



e^^p^k^K^'qPQ'^ = -XMBXK^q^m^, sin sin sin , 

01 {6k) is the angle between the i~ {K) three-momentum vector in the i'^i" (Kit) rest frame 
and the direction of total {K*) three-momentum vector defined in the B rest frame, cp 
is the angle between the normals to the planes defined by i'^i' and the Kn, in the B rest 
frame. The differential decay rate is then given by 

dT = \ f \M\^XXKXedq'^dcos9Kdcos9id<p , (10) 

assuming a narrow width approximation for the decay K* — > Ktc. 

It can easily be seen from eq.@ that, the only terms proportional to sin((y9) or sin(2y9) 
are those that depend on the imaginary parts of the products any two of a, (3 or p. For 
instance only the coefficient of Im(ai/32 -|- ctRPji) is proportional to sm{2ip). Hence we can 



7 



isolate this term by considering the following asymmetric width in terms of the differential 
decay rates of the B meson with respect to ip, 



- "2^. dT 



The imaginary part in the term under consideration can be due to either a strong phase 
or a weak phase. Such CP violating asymmetries can be obtained not by considering the 
difference of differential rates for B and B , but the sum of these rates. It follows trivially 
from eqn.(|^ and ^ that the asymmetric width for B{B ) is, 

Ari(Af oc ±E{l«ill&Ll sm((5f ) ± (01^)) + L^R}, (12) 



where = {62' — 5^) and (j)-'^ = (0^^ — S^''). We define the asymmetry Ai as the sum of 
the asym 
we have, 



All + All „ 
the asymmetric widths, normalized to the total widths, Ai = — - — ^ — and from eq.(|l^), 



Ai oc cos(<5f ) sin(0i^) + L R} , (13) 

which is nonzero if and only if there is CP violation represented by non-zero phases 
|T^,2D.2^|. 6 can arise from electromagnetic final state interactions, which are negligible 
and ignored. For top quark in the penguin loop, the case of 6 being nonzero due to inter- 
mediate quark on shell, does not arise. 

It is also possible to construct a different asymmetry that isolates another combination 
of the imaginary terms. Such an asymmetry considers the difference distribution of the 
same hemisphere and opposite hemisphere events, and the asymmetric width in this case 
can be defined by, 

AT2 = { - )d^ dcosOi / c/cos^^- — — - , (14) 

JO Jn Jd Jd d cost! id cos UKd(p 

f AFa + AT2 

where / = / — . Analogous to Ai, we define the asymmetry A2, A2 = = — . 

Jd J-1 Jo T + T 

The asymmetries Ai and A2 are evaluated to be, 

Ai = -2xA / dq'^C{aoV - Aog)X\ A2 = x f dq'^CF ^-^X^ , (15) 

J J niK* V g 



where, 

_ a^G% nib 

^ ^ 24 tt^Mb (r + f ) (Mb + tuk* ) ' 
C = {kC'^ + \C^)Cl^, A = Ml -ml., 

F = {2X''Ml{a+ V - A+g) + {aoV - Aog) Ak ■ q} . 

Note that these asymmetries are independent of Cg^ . Since, Cg''** terms are suppressed 
m comparison with C'^'^, for all our estimations, we ignore them. 

For a rate asymmetry A to provide a n standard deviation signal of CP violation, we 
require that A = n/y/N, where N is the total number of events in the channel. Thus for 
A2, at la (1.64 a) level, it is possible to place the bounds, < 0.34 (0.92) for A = and 
|A| < 1.02 (2.75) for k = 0, using 10,000 [B K*e+e-) events. This corresponds to 2x10^ 
i^'s, however, by adding the {B K*^'^^~) mode, 10^ i^'s will be required. Much weaker 
bounds would be obtained from Ai. 

However, we can improve the statistical significance, by looking at the distributions, 
shown in Fig. 1. We bin the asymmetry A2 (Ai being small is ignored, and we shall drop 
the subscript 2 hereafter) and use a fit to obtain bounds on the parameters k and A. The 
asymmetry as a function of is given by, 

(/ -/ )d(p d cos 9i I d cos 9k Tsumiq^) 

Mq') - Z X — # > (16) 



/■TT /'27r /■! /■! ' 

(J^ + J J ^dcos9i J ^dcos9KVsum{(f') 



where Vsumiff) = -m 7n n — '^(^ + T) is the sum of the differential widths of B 

dq^ d cos 9 1 d cos 9Kdip 

and B. The average value of asymmetry in the i*'^ bin is then, 

I dq^A{q^) \ ' , (17) 



' Tb Ni Jq^ . ^ dq'^ 

where, Nb is the total number of S's, is the total B width and q^in,Qmax minimum 
and maximum q'^ values in the bin; the number of events in the i^^ bin are, 

Tb Jq" . dq^ ^ ' 
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The value of the asymmetry, in the i*'* (f bin, coming from k and A contributions, is of 
the form 

A = (19) 

In Xj, terms quadratic in k and A are ignored, which is a reasonable approximation as long 
as A are small. The observed asymmetry Af^ in the i^^ bin is assumed to be chosen from 
a Gaussian distribution with mean equal to the theoretical asymmetry A*'' and variance af. 
Then, the method of least squares gives, 

nhin / Aobs AthX"^ 

x^ = EF^-r^ • (20) 

i=l \ J 



The statistical errors are taken to be CTj = yl/N. We evaluate the difference, Ax^ = 
XsM ~ Xmim where Xsm ^^e corresponding to the SM values of the parameters i.e. 
(0,0) and xLin evaluated from the values {k,X) that minimize x^. The linear dependence 
on the parameters k, A results in the simple form. 



2 / ~ \ 2 

/ utti \ ( Xbi \ 2\kaih. 



- 



Ax' = E ^ ^ ^ ■ (21) 



i=i '^i 



Now, Ax' — r? will give the n-standard deviation bounds for k and A. 

For 20 bins, for 10,000 {B — > K*e^e'^ events, it is possible to obtain the improved 
individual bounds, 

\k\ < 0.25(0.42) at 68.3(90)% C.L. 

|A| < 0.76(1.25) at 68.3(90)% C.L. (22) 

The possible combined bounds on k and A are shown in Fig. 2. 

To conclude, we have studied CP violating anomalous couplings k and A contributing 
to the process B — > K*t^l~ — > K7ri'^i~ . These anomalous couplings can be constrained by 
constructing an asymmetry, requiring the addition of B and B events. No strong phases, 
nor any flavour/time tagging are required for this technique. Although, EDM of neutron 
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places strong bounds on the couplings R and A, the CP odd asymmetries studied here, 
depend on a different combination of these anomalous couplings and thus provide useful 
additional information. A being the coefficient of a higher dimensional operator, is expected 
to be better constrained by high energy collider events than from rare B decays. The recent 
results of DO provide good bounds on A, but for k the approach discussed here can provide 
a much tighter constraint. 

We would like to thank G.Date and G.Rajasekaran for discussions. 
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FIGURES 

FIG. 1. (a) The differential branching fraction as a function of . (b) The asymmetry ^2(9^) 
as a function of c^. The sohd line is for k = 1(A = 0) and the dashed for A = 1(k = 0). 

FIG. 2. Limits on the CP -violating anomalous coupling parameters h and A. The region 
inside the solid (dashed) lines represents the 68.3(90)% C.L. limits. 
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TABLES 

TABLE L Relations between the form factors used in this paper, a quark model (QM) that 

k ■ K 

reproduces heavy quark limit . = {K^ — C^/x), o"^ = 967r^/(m|^, Ai^-'^), Q = — ^ — i and A 



are defined in the text. 



QM i 



m 



K* 



A -25cr 
B ao A cr 

C 2a+W -qa -CB 



V -2 ^ a 

Mb + rriK* 

Mb + rriK* 

W-qa-C£ 



Mb + tuk 



15 




Figl 




Fig 2 



